Kondo-transport spectroscopy of single molecule magnets 
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We demonstrate that in a single molecule magnet (SMM) strongly coupled to electrodes the Kondo 
effect involves all magnetic excitations. This Kondo effect is induced by the quantum tunneling of 
the magnetic moment (QTM). Importantly, the Kondo temperature Tk can be much larger than the 
magnetic splittings. We find a strong modulation of the Kondo effect as function of the transverse 
anisotropy parameter or a longitudinal magnetic field. For both integer and half-integer spin this 
can be used for an accurate transport spectroscopy of the magnetic states in low magnetic fields 
on the order of the easy-axis anisotropy parameter. We set up a relationship between the Kondo 
effects for successive integer and half-integer spins. 

PACS numbers: 72.10.Fk, 75.10.Jm, 75.30.Gw, 75.60.Jk 



Introduction. — Single molecule magnets (SMMs) al- 
low the study of quantum phenomena on a mesoscopic 
scale, namely the quantum tunneling of the magnetic mo- 
ment. SMMs such as Mni2 or Fes have attracted intense 
experimental and theoretical investigation in recent 
years. Due to weak intermolecular interaction molecu- 
lar crystal properties can be assigned to single molecules 
described by a large spin (S > 1/2), and easy-axis and 
transverse anisotropics. Recently, two groups 0,131 have 
trapped a single molecule magnet (Mni2) in a three ter- 
minal transport setup and measured transport through 
the single molecule. Electron transport fingerprints due 
to both sequential 0, Q and inelastic co-tunneling Q 
processes were observed and associated with the molecu- 
lar magnetic states. Theoretical works [JIB-El predicted 
that fingerprints of magnetic quantum tunneling (QTM) 
can in principle be identified in transport measurements 
in the charge (sequential tunneling) as well as in the 
spin fluctuation (Kondo) regime. Especially the latter 
regime of strong coupling to the electrodes is of impor- 
tance since the Kondo effect shows up as a sharp zero-bias 
anomaly with width given by the Kondo temperature 
Tjc ■ The Kondo effect has been observed experimentally 
in many other systems with small spins ( e.g. quantum 
dots 0E HE EH and sin S le molecules USUI). For 
SMMs with half-integer spin S in zero magnetic field it 
was shown [j| that the Kondo effect arises from a cooper- 
ation of both spin exchange processes with the reservoirs 
and the intrinsic tunneling of the magnetic moment. The 
strong coupling fixed point is of pseudo-spin- 1/2 type. 
In this limit the exchange coupling J is weak, such that 
the resulting Kondo temperature Tj{(J,A) <C A, where 
A < 1 meV is the scale of the anisotropy splittings be- 
tween the magnetic states. In this case, the Kondo effect 
involves only the two magnetic ground states of the SMM 
and does not occur for integer spin. 
In this Letter we study the Kondo effect in the ex- 
perimentally more favorable regime of strong exchange 
(J ~ 0.1 eV) where the Kondo temperature becomes 
larger than the magnetic splittings Tr-(J, A) 3> A. We 
find that excited magnetic levels on the SMM, belonging 
to different topological sectors (with respect to rotations 



around the z-axis), become essentially involved in the 
Kondo effect. In zero magnetic field and for half- integer 
spin the Kondo effect can be modulated by changing the 
transverse anisotropy, resulting in a sequence of Kondo 
effects associated with different magnetic excited states. 
A longitudinal magnetic field can induce two important 
new effects: (1) For both integer and half-integer spin 
the Kondo effect is suppressed at each anticrossing of 
magnetic levels belonging to the same topological sector. 
The corresponding sharp magnetic field scale is deter- 
mined by the transverse magnetic anisotropy. Therefore 
the study of the magnetic field dependence allows for 
an accurate spectroscopy of the magnetic states. This 
modulation of the Kondo effect allows for an experimen- 
tal proof of the existence of QTM in a single molecule 
in a transport junction and to determine the important 
microscopic parameters characterizing the SMM in that 
setup. The magnetic field induced Kondo effect for inte- 
ger spin SMMs allows many molecular magnets like e.g. 
Fes to be studied without charging the molecule to obtain 
a half-integer spin. (2) We find that the Kondo effects for 
successive integer and half-integer spin S display a close 
correspondence when shifted in magnetic field energy by 
the easy-axis anisotropy parameter D. 
Model. — We consider a SMM (Fig. ^| in a transport 
setup where the applied voltages, charging effects and low 
temperature suppress single-electron processes 5j . In the 
presence of a longitudinal magnetic field the Hamiltonian 
reads H = H S mm + H ex + i? res : 



-ffSMM 
Hex 
Hrcs 



-DS 2 Z 
JSs 



\B 2 (S\ + S 2 _) + H Z S Z 



E 



e kaa kcr aka 



(1) 

(2) 
(3) 



where S z is the projection of the SMM's spin on the 
easy z-axis, and S± = S x ± iS y . The terms in Eq. 
describe, respectively, the easy-axis magnetic anisotropy 
of the molecule, the transverse anisotropy perturbation 
and the coupling to a magnetic field (H z ) along the easy 
axis. The g-factors are absorbed into the magnetic field. 
For simplicity we have taken a 2-fold rotation-symmetry 
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FIG. 1: Energy diagram for Hm with B2 <C D without (left) 
and with (right) magnetic field. The two subsets of eigen- 
states of different magnetic symmetry are indicated by full 
and dashed bars. 

axis, which is dominant in many molecular magnets. The 
Hilbert space of -ffsMM is split into two disjoint topolog- 
ical subspaces a = ± since the 2-fold rotational symme- 
try about the z-axis is preserved under the longitudinal 
magnetic field. While for half-integer S these spaces are 
spanned by an equal number of ra+ = n_ = 5+1/2 
basis-states {\S,M)} M=<7S ,a(s-2),...,-<7(S-3),-a(s-i), for 
integer spin they are spanned by n+ — S + 1 states 
{\S, M)}m=-S,-S+2,...,S-2,S for a = + and by n_ = S 
states {\S, M)} M =s+i ,-s+a,...,s-a,s-i for o = -. The 
eigenstates \lo) of H^mm with energies E\ a are labeled 
by I = 1, 2, ..,n CT in order of decreasing energy in each 
subspace (see Fig. Q). Eq. J5J) describes the exchange 
coupling of the molecular spin to the effective reservoir 
Eq. (bandwidth 2W and constant density of states 
p). The electronic states are labeled by ko and denote 
the even combination of left and right physical electrode 
states |3 113 • The local electron spin in the reservoir 

is s = J^kk' J2aa' a \a T <y<y' a k'a' /2 where r is the Pauli- 
matrix vector. The exchange coupling J is induced by 
virtual electron tunneling processes and is antiferromag- 
netic due to the strong energy- and charge quantization 
effects. The Kondo spin-scattering of conduction elec- 
trons off the SMM transfers charge between physical left 
and right electrodes and can resonantly enhance the lin- 
ear conductance. The half-width at half-maximum of the 
resulting zero-bias differential conductance peak at T = 
is the Kondo temperature Tr. 

For a qualitative discussion it is useful to change to the 
exact representation of Eq. (|2J) in the eigenbasis of the 
molecular states \lu): 

H «* = E ( E J «' P »'^ + E J n'« 1 s * I ■ ( 4 ) 

\i—x,y er / 

The Hamiltonian Eq. (J2J describes the exchange pro- 
cesses involving reservoir electrons and internal mag- 
netic degrees of freedom of the SMM. This projection 
makes explicit that the QTM parameter E>2 and the 
external magnetic field H z modulate the effective cou- 
plings through the matrix elements of the molecular 
spin operator S. The first, transverse, term describes 
the spin-scattering involving a pair of states \l+) and 



\V— ) in terms of pseudo-spin-1/2 operators Pi Vl iPf v = 
- I ± \l'-){l + |)/2 with effective exchange cou- 
plings 

J^{ y = J{l + \S+±S-\l'-). (5) 

In total, there are n_|_n_ such pairs of states from op- 
posite topological sectors. The longitudinal couplings 
read Jf Va = J (la\S Z \V a) . Longitudinal spin operators 
Pf v = + I - - |)/2 can only be introduced 

in a unique way if an approximative projection onto a 
single pair is made. For instance, one recovers the pro- 
jection of Ref. Q onto the ground state pair for zero-field 
and half-integer S by using time-reversal symmetry and 
by truncating the excited states. However, such a trun- 
cation is not valid in the regime of interest here: the 
strong exchange coupling J gives rise to a Kondo tem- 
perature that can be larger than the magnetic splittings 
(T K > 2SD for B 2 < D). The Kondo effect can occur 
irrespective of whether the pair of states \l+) and \V— ) 
are ground- or excited states of the SMM and whether 
they are exactly degenerate or not. The Kondo effect 
therefore involves contributions from multiple topologi- 
cal pairs. This we confirmed explicitly by calculating the 
projection of the full many-body ground state onto each 
molecular eigenstate using the NRG. Thus the full ex- 
pression for H CXl Eq.(0J), must be retained for which a 
unique decomposition into a sum of independent pairs is 
not possible. This is due to the fact that a single eigen- 
state of the SMM \la) is paired with n_ CT > 2 states 

\V-v)- 

Physical picture. — Variation of either the QTM parame- 
ter E>2 or the magnetic field H z results in anticrossings of 
magnetic levels. As shown in Fig. ^ the variation of the 
magnetic field leads to an alternating sequence of degen- 
eracies between pairs of levels which are in different sub- 
spaces (crossings) or in the same subspace (anticrossing) . 
Importantly, this occurs in the low magnetic field en- 
ergy window where the Zeeman splittings are still smaller 
than Tk {Ei a — Eir a > < Tk)- The effective couplings in 
Eq. are modulated strongly at each anticrossing lead- 
ing to a suppression of the Kondo effect on a scale ~ E>2 , 
as we now explain. At an anticrossing states \l<r) and 
\l'a) from the same topological sector are close in energy 
and strongly hybridize due to the transverse anisotropy. 
Following two levels adiabatically during an anticrossing 
they interchange their role and one basis state picks up a 
relative phase tt. For example, in Fig. ^thc level |1+) is 
moving upwards and the level |2+) downwards in energy 
after the anticrossing. For each topological pair involving 
one of the anticrossing levels this leads to a sign change 
of one of the transverse couplings, as can be seen easily 
from Eq. (J5J . The latter therefore vanish at the anticross- 
ing and the Kondo effect is expected to be suppressed. 
This happens on a magnetic field scale proportional to 
the tunnel splitting, i.e. will occur as a sharp feature 
in the linear conductance as function of a longitudinal 
magnetic field. In a similar way, in zero field the varia- 
tion of the QTM parameter P>2 itself leads to a series of 
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anticrossings, which for a half-integer spin preserves the 
2-fold Kramers degeneracy of all levels. This happens 
due to the non-uniform level spacing and the coupling of 
states occurs in different orders of perturbation theory in 
£?2 (higher-lying levels hybridize more strongly for weak 
B 2 ). 

Method. — We use Wilsons's NRG [H El to treat non- 
perturbatively the full model H = Hqmm + H cx + H res . 
As NRG input parameters we use number of states 
N = 1000, discretization A = 2 and D = 5 x IQ- 5 W . 
We analyzed the NRG level flow as function of it- 
eration number iViter in order to determine the low- 
temperature fixed point for H z = 0. The Kondo temper- 
ature is defined as the energy scale where the crossover 
to strong coupling takes place. We also calculated the 
spectral function within the T- matrix approach |18| : 
A a (uj) = — —ImT a (uj + iS) where the T-matrix is 
T a {u) = {{oJ,Ol)) with O a = !(c ,_ CT S- CT + <rco,*S M ) 
and co,ct = X)fe c k.& is the electron operator on the first 
site of the Wilson chain. The low temperature conduc- 
tance G = —{e 2 /2)J2af duA a (uj)df(tu)/dtu is propor- 
tional to the spectral function, where f(u>) — l/(e u ^ T + 
1). The effects considered here are difficult to capture in 
a poor-man's scaling approach due to the strong coupling 
and the many excited states involved. However, another 
useful guide to the full NRG results in this regime is 
the spin binding energy AE, obtained by diagonalizing 
exactly the molecular Hamiltonian Eq. coupled by 
exchange to a single conduction electron spin. This is 
similar to the first step of a NRG calculation and will be 
denoted as the zero-bandwidth model. The ground- to 
excited-state gap AE thus obtained follows the modula- 
tion of Tk by £?2 and H z accurately, although the scales 
of AE and Tk strongly differ. Apparently, for a siz- 
able range of the relevant exchange strengths J < W 
the zero-bandwidth estimate AE is renormalized uni- 
formly by the coupling to the remaining conduction band 
electrons. This estimate breaks down, when the modu- 
lation involves a suppression of Tk far below the scale 
of the magnetic splittings: then the renormalization be- 
comes non-uniform and the variation may deviate from 
the NRG. 

Kondo effect due to excited states tuned by QTM. — We 
first consider the case of zero magnetic field. For integer 
spin this corresponds to an anticrossing: the eigenstates 
are paired to nearly degenerate states from the same sub- 
space which arc splitted by the transverse anisotropy. 
Consequently, as discussed above the Kondo effect is 
suppressed, resulting even in a dip in the spectral func- 
tion, see below. In contrast, for half-integer spin many 
topological pairs are degenerate (crossing) due to time- 
reversal symmetry, see Fig.^ The NRG converges to the 
spin-1/2 strong coupling fixed point, indicating a com- 
plete screenin g of t he magnetic degrees of freedom. As 
shown in Fig. |2(a)| the Kondo temperature shows an os- 
cillatory dependence on B2/D. The number of oscilla- 
tions for < B2 < D is S — 1/2. When decreasing 
the exchange coupling J the smallest peaks disappear 




FIG. 2: Dependence of Tk on the transverse anisotropy 
B2/D: (a) Tk from the full NRG calculation, (b) spin binding 
energy estimate AE. Parameters: S = 7/2, D = 5 * 10~ 5 W. 



first, leaving only the monotonic increase of Tk up to the 
broad maximum centered at B2 ~ D in the limit where 
the magnetic excitations can be neglected 5J. The spin 
binding energy AE captures this dependence on B2/D, 
2(b)| Interestingly, we find that the 



as shown in Fig. 

nature of the ground state of the zero-bandwidth model 
changes with each oscillation which can in turn be re- 
lated to an anticrossing of magnetic states on the SMM. 
For ins tance , for weak B2 <C D near the lowest peak 
in Fig. 2(b) the highest excited Kramers doublets are 



strongly mixed into the zero-bandwidth model ground 
state. In contrast, for B2 ~ D near the highest peak in 



Fig. 2(b) the ground doublet dominates in the ground 
state. This demonstrates that the observed strong cou- 
pling Kondo fixed point in the NRG indeed originates 
from a screening of magnetic degrees of freedom involv- 
ing excited magnetic states of the SMM. Which of these 
excited states are important depends on B 2 /D. 
Kondo- spectroscopy of SMMs. — We now focus on the 
most relevant case of fixed weak QTM (B2 <C D) and 
vary the longitudinal magnetic field H z . For many mag- 
netic field values the spectral density shown in Figs. 3^ 
3(d) displays a zero-bias Kondo resonance. The corre- 



sponding many-body ground state found in the NRG 
at the strong coupling fixed point is non-degenerate. 
This peak is strongly modulated at avoided crossings of 
states from the same subspace at H z /D w 2k (2k + 1), 
k = 1,2,... for integer (half- integer) spin. As discussed 
above, this is related to the strong suppression of effective 
exchange couplings on the magnetic field scale B2 near 
an anticrossing. In contrast, the spectral function varies 
smoothly at level crossings located at H z /D sa 2k + 1 
(2k) for integer (half- integer). In this way the two scales 
D and B2 can be identified in the magnetic field depen- 
dence of the linear conductance. For integer spin S = 1 
the spectral function in Fig. 3(a) shows a dip for H z = 
(anticrossing of |1+) and |2+)). At fields H z ~ B 2 a 
Kondo peak is induced which reaches maximal height at 
H z = D (crossing of |1+) and |1— )). This Kondo peak 
is subsequently suppressed and splitted due to the Zee- 
man level shifts, see Fig.^ For half-intege r spin 5 = 3/2 
the zero-frequency peak at H z = in Fig. 3(b) is broad- 
ened with increasing H z . This broadening is maximal at 
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FIG. 3: (Color online) Spectral function, qualitatively equiv- 
alent to dl /dVbias as function of H z /D and frequency ui (a) 
for smallest non-trivial integer spin S = 1, (b) half- integer 
spin S = 3/2, and (c) large integer 5 = 3, (d) half-integer 
spin S = 7/2. Parameters: J = 0.15W, D = 5 * 10 -5 W and 
B 2 /D = 0.1. 



H z — D and is also captured by the H z /D dependence of 
the spin-binding energy AE (not shown). The full NRG 
result in addition shows a sharp dip superimposed on the 
broad peak at H z = D (anticrossing of |1+) and |2+)). 
At this point the situation is very similar to the S = 1 
case, provided one neglects the contribution of the low 
lying state 2— ) and reduces the magnetic field by D, see 
Fig. n Indeed, increasing the field further to H z = 2D, 
the Kondo effect reappears precisely as for S = 1 near 
H z = D. For larger spin the physics is qualitatively 
the same: the spectral function in Figs. 3(c) 3(d) dis- 



plays several dips of the Kondo peak close to each anti- 
crossing. The correspondence pointed out above is more 
general: the linear conductance of a SMM with spin S at 
magnetic field H z corresponds qualitatively to that of a 
SMM with reduced spin S — 1/2 and reduced magnetic 
field H z — D. This can be shown using the Hamilto- 
nian Eq. 10} in the molecular eigenbasis. For fixed J the 
correspondence in the conductance is most obvious for 
subsequent S values. The conductance for spin S near 
H z = nD (n = 1, 2, . . . , 25 — 1) can even be compared 
with that for S — n/2 around zero- field by iterating n 
times, which we have checked. However, a clear corre- 
spondence appears only when J is properly adjusted. 
Conclusion. — Magnetic field Kondo transport spec- 
troscopy can determine the absolute value of the spin 
and the magnetic parameters of a SMM in a transport 
junction. The strong Kondo effect in SMMs discussed in 
this work can be accessed experimentally by increasing 
the exchange coupling. In STM setups one can reduce 
the distance to the molecule, or change the molecular ge- 
ometry by a bias- voltage pulse |l9j . In 3-terminal mea- 
surements 0, El one can tune the exchange coupling with 
the gate voltage: J <x 1/\V* — V g \. In addition, both the 
charge- and spin- state can be changed by tuning the 
gate voltage to opposite sides of the charge degenerac 
point V g = V* . If a single-electron transport current 
is observed at this point c.f. 2, 3], S only changes by 1/2. 
Then the predicted clear correspondence between the lin- 
ear conductance for subsequent integer and half-integer 
spin values provides an additional check on the physics. 
Importantly, the spectroscopy can be done at tempera- 
tures above the energy scale of the magnetic splittings 
and requires only low magnetic fields which do not de- 
stroy the Kondo effect by Zeeman splittings. Finally, we 
have checked that small fixed transverse magnetic field 
perturbation does not destroy the reentrant behavior of 
the Kondo effect as the longitudinal field is varied. 
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